Vortex dynamics in two-dimensional Josephson junction arrays 
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The dynamic response of unfrustrated two-dimensional Josephson junction arrays close to, but 
above the Kosterlitz-Thouless(ifr) transition temperature is described in terms of the vortex di- 
electric function e(aj). The latter is calculated by considering separately the contribution of ^CfreeS> 
vortices interacting by a screened Coulomb potential, and the <Cpair motions of vortices that are 
closer to each other than the KT correlation length. This procedure allows to understand various 
anomalous features in e(o>) and in the flux noise spectra that have been observed experimentally 
and in dynamic simulations. 
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Two-dimensional (2c?) Josephson junction arrays 
(J J A) are described by the classical XY- model in which 
the superconducting phases Oi of neighboring sites I are 
coupled by the Josephson interaction. The relevant col- 
lective phase excitations are vortices (V) and antivor- 
tices (A). They behave as a 2d neutral Coulomb gas 
(CG). Their charge is related to the Josephson coupling 
J by q 2 Q = 2ttJ . At the Berezinsky-Kosterlitz-Thouless 
(BKT) transition the CG goes over from a <C dielectrics 
phase consisting of T^A-pairs to a <CmetallicS phase con- 
taining free V and A. One of its main signatures is 
the universal jump of the helicity modulus at Tbkt- It 
has been observed experimentally, both in the current- 
voltage characteristics of the array and in its response to 
an electromagnetic field of very low frequency [1,2]. 



The dynamics of J. J As is less well understood. For 
the equations of motion for the superconducting phases 
the -Cresistively shunted junction models (RSJ) can be 
used for arrays in which electrostatic charging effects may 
be neglected. Based on this model an equation of motion 
for the vortex excitations can be derived [3], describing 
them as massless point particles, subject to a friction 
force coming from normal current losses and interacting 
via the 2d Coulomb interaction, varying essentially with 
the logarithm of their distance. The dynamics of such 
a system is usually treated by considering separately the 
motion of bound pairs and of free particles. The simplest 
picture for free motion is given by the Drude form of the 
dynamic dielectric function e(uj), expressed in terms of 
the friction constant. Minnhagen [4] has developed a 
more sophisticated expression for e(u>) by assuming that 
this quantity can be derived from the static wave num- 
ber dependent dielectric function, taken to have a Debye 
screening form, by replacing the wave number by fre- 
quency. The result - usually referred to as <cMinnhagen 
phenomenology S (MP) - differs in an essential way from 
Drude's (D) behaviour. In particular Re(l/ cd(u)) oc uj 2 , 
whereas Re(l/eMp(^>)) oc \ui\, and the so-called <Cpeak- 



ratioS r, given by 



where Umax is the frequency at which Im(l/e) has its 
maximum, is rjj — 1, whereas tmp = 2/7T is smaller. 

The dynamics below Tbkt is usually treated by aver- 
aging the dynamic response of a pair of separation d over 
a probability distribution for d [5,6]. 

Three main types of experiments aim at elucidating 
dynamical properties of J J As. The exponent of the 
non-linear current-voltage characteristics is related to the 
dynamic critical exponent describing the critical slow- 
ing down of the 2d CG near the BKT transition (see 
[1]). The dynamic conductance G(u>) of the array can 
be inferred from measuring the dynamic response of the 
array to a time dependent current in a two-coil experi- 
ment [1,2,7]. Measurements on frustrated arrays [8] have 
produced results that are closer to the MP prediction 
than to the simple Drude form, in particular as far as 
the above-mentioned peak-ratio and the frequency de- 
pendence of l/e(u>) are concerned: Re(l/e(u>)) cx \lo\ over 
a sizable range of frequencies, which is a signature of 
anomalous MP dynamics. Similar results have been ob- 
tained by analytical calculations [6,9,10]. 

Measuring the temporal fluctuations of the magnetic 
flux through a given area of a J J A [11,12] gives another 
insight into the dynamics of the currents in the array, and 
thus of the vortex system. At sufficiently low frequen- 
cies, above Tbkt, the time Fourier transform S c / > (oj)of 
the noise is white (frequency independent) . For larger uj 
these experiments show a rather extended region where 
S<j,(us) oc 1/u), which is unexpected and needs a deeper 
explanation. 

Various analytical and numerical methods have been 
used aiming at explaining the somewhat unexpected 
<C anomalous S dynamics of J J As revealed by the above 
mentioned experiments. Dynamic simulations for J J As 
[13-16] have been performed based on the equations of 
motion for the phases of the array or for the CG [17]. 
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The extensive work of Minnhagen and collaborators has 
confirmed the existence of a frequency interval in which 
the anomalous MP dynamics can be seen, both for zero 
field and for frustrated arrays. On the other hand the 
flux noise spectrum obtained in these and other calcu- 
lations usually only show a common l/u> tangent to the 
curves for different temperatures, but no extended region 
with l/w-noise. For higher frequencies S^,(uj) cx 1/w 3 / 2 
[18-20], which is characteristic of vortex diffusion or 
S<f,(oj) oc 1/w 2 , for even higher frequencies [20]. Other 
calculations do yield l/w-noise [21] or cx \/uj a ^ 

with an exponent a(T) close to 1, but slowly varying with 
temperature [22]. 

Our theoretical approach to J J A dynamics is based on 
the equations of motion for the 2d neutral CG containing 
a friction term and the effect of the Coulomb interaction. 
The key quantity for calculating all observables of inter- 
est is the dynamic charge correlator: 



>(k, 



dte~ zt (p(k,t)p*(k)) 



(2) 



p(k) being the Fourier transform of the charge density. 
The equal time correlation function S(k) — { | /o(k) | 2 ) is 
the charge structure factor of the CG. From 4> pp one can 
obtain the dynamic charge susceptibility and the dielec- 
tric function of the CG [25] . The later can be expressed 
in terms of a dynamic vortex mobility // 



e(w) = 1 - 



2tt<7oM(w) 



(3) 



and it yields the superconducting part of the dynamic 
conductance of the J J A: 



iu>e(u) 



(4) 
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FIG. 1: Re[l/e(u)] and -Im[l/e(u))] vs u/u a plot (dot- 
ted curves for —Im) for free V and A. From more 
thickness to less for T/Tbkt— 1.01, 1.05, 1.10. 

Assuming that supercurrent is the main source of the 
flux threading through a coil of radius R at a distance d 



from the array plane, S^(u>) is related to 



2(2^ 2 i?/i ) 2 (^) 2 (JA) 2 x 



2„-2kd 



dk 



Ji(kR) 2 e 
k(l + Xk) 2 



-Recf) pp (k, -iw) (5) 



Ji being the first order Bessel function and A the mag- 
netic penetration depth of the J J A. 

Our goal consists in elucidating - by analytical means 
- the critical dynamics of the I^A-system when Tbkt 
is approached from above. In this regime two types of 
motion, occurring at different length scales, have to be 
considered simultaneously. Particles closer together than 
a typical T-dependent critical distance d c are supposed 
to move as a bound pair, at least for times smaller than 
some typical life time r (depending also on T). They 
interact through the Coulomb potential, screened as in 
a dielectric. On the other hand, particles further apart 
from each other than d c will be supposed to move like 
^Cfree» - i.e. <Cunbound» - particles, subject to a 
metallically screened potential. A complete theoretical 
approach should, of course, treat the two types of time 
evolution with the same mathematical tools. Here we 
have chosen to evaluate separately the mobilities /if, for 
<Cbound3> and p,f for <g;free;§> motion. The two results 
are then combined by calculating the total mobility as a 
weighted sum of Hb and [if , the relative weights being 
given by the density of paired and free excitations. 

For the free motion /i/(w) is found by applying Mori's 
technique of calculating dynamic correlation functions to 
4> pp [6,10]. As usual higher order correlators, showing up 
when the equations of motion for the individual particles 
are used, are factorized. In order that this factorization 
be adequate for a system with long range interaction, we 
replace the Coulomb potential by its screened version. 
Its Fourier transform has the form 



V 3C (P) 



to 



r 
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(6) 



with £ sc being the (metallic) screening length which we 
identify with KT correlation length. This replacement is 
standard in the treatment of quantum Coulomb systems 
[23]. The resulting density correlator (2) has the usual 
form 



• J i>p 



CM) 



S(k) 



z + 



k B Tk 2 
S(k) 7 (z) 



(7) 



The generalized friction function 7(2) = /i _1 (z) is the 
memory kernel of the charge current density. 'Jf(z) for 
free motion shows MP behavior over a frequency range 
which increases when Tbkt is approached. There are 
three frequency regimes, separated by the scale frequen- 
cies uj a — with a being the lattice constant and 

w £ = rf^"' -^ e [T/( — ^ UJ )] i s fl a t f° r w < w £> ^ coincides 
with the bare friction parameter T for lu > uj a , whereas 
in between i?e(7/) oc In \oj\. This is precisely what MP 
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predicts, and the corresponding inverse dielectric func- 
tion shown in figure 1 indeed follows the MP prediction 
for a frequency range that increases when Tbkt is ap- 
proached. The values of the peak ratio (1), varying be- 
tween 0.67 and 0.73 for the three temperatures shown 
in figure 1, are also much closer to MP than to Drudc. 
Thus, in this framework the anomalous MP dynamics is 
explained by the increasing influence of the long range 
Coulomb force, which makes the motion more and more 
<Csluggish;^. Whereas at high temperatures the poten- 
tial is (metallically) screened for all relevant length scales, 
this screening becomes less and less efficient when Tbkt 
is approached where the screening length £ sc diverges. 
However, contrary to MP the response always crosses 
over to Drude-like when uj < lu^ . If this were not the case 
the arrays would still be superconducting above Tbkt 
[10], whereas in reality a finite flux-flow resistance yields 
a finite conductance G(u> = 0). The corresponding flux 
noise spectrum, however, does not show any extended 
1/ui region. Both results, l/e(u>) and S^ui), are very sim- 
ilar to Minnhagen's simulations. One may thus draw the 
conclusion that anomalous MP behaviour for the dielec- 
tric function does not yield the - equally <C anomalous 3> 
-1/uj noise, which should thus have a different origin. 
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FIG. 2: S<f,(ui) vs u>/uja plot for V A pair dynamics in 
method (b) with fixed fc ol = 1.0. From more thickness 
to less for T/T B kt= 1.01,1.05,1.10. 




FIG. 3: S,/, (oj) vs iojuja plot with T dependent fc„2 for V A 
pair dynamics in method (a). From more thickness 
to less for T/Tbkt= 1.01,1.05,1.10. 

Turning now to the <Cpaired^> motion we average the 



dynamic polarizability of a single pair over a suitable 
probability distribution function 4>(r) for the pair size r 
[5,6]: 



rdr- 



(r) 



We have chosen a normalized Boltzmann factor for <h(r 



(8) 



r) oce p " e m 9(d c - r) 



(9) 



using two different forms for the static dielectric func- 
tion screening the bare Coulomb potential ; (a) e(r) fixed 
at e(oo, Tbkt) = /{/HzbTbkt) and (b) e(r) length 
scale dependent, as determined in the framedwork of the 
KT scaling equations [4]. The distribution is cut off at 
the length scale d c . For (a) this is the BKT correla- 
tion length £;bkt, and for (b) it is the distance where 
the screened potential turns over from attractive to re- 
pulsive. Both have the characteristic BKT temperature 
dependence with simply somewhat different parameters. 
The two procedures (a) and (b) yield similar results. In 
the same spirit the force constant in the denominator of 
(8) is given by 



m = 



to 



2 e(r) 



(10) 



We have also introduced a finite pair life time [6] given 
by 



1 



r(r) 



j exp{ 



ql \n{d c /a) q 2 Q ln(r/a) 



e(d c ) 



e(r) 



} 



(11) 



where u> , a free parameter, is an <C attempt frequency 3> 
for pairs trying to escape the barrier up to the maximum 
of the potential at the distance d c by thermal excita- 
tion. The inverse dielectric function resulting from (8) 
shows again marked deviations from Drude behaviour, 
in particular for i?e(l/e(w)) varying like w s ( T ) for inter- 
mediate uj. The T-dependent exponent s(T) « 1/6 to 
1/3 is due to the r-integral (8) where the dependence of 
4>, f and 1/r on r imposes a particular exponent to the 
frequency dependence of the susceptibility. This even 
<Cmore anomalous^ frequency dependence is reflected 
in the peak ratio r (equation (1)) which is now even 
smaller than the MP value of 2/tt. The low oj cross- 
over to Drude behaviour occurs at a critical frequency 
u) c , which is found to, vary as d~ 2 , whereas for uj > uj a 
it crosses over to the high frequency Drude form. In the 
window uj c < uj < ui a the behaviour is anomalous. This 
has very interesting consequences for the flux noise spec- 
trum presented in figures 2 and 3. The values used for 
the parameters showing up in the integral (8) are given 
in the figure caption. The structure factor in the density 
correlator (7) has been given the following form 



S(k) 



(12) 
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which respects the charge neutrality of the V — A sys- 
tem and the correct limit for large k. Two different ideas 
have been used for the inverse characteristic length k a 
in S(k). For 2d neutral Coulomb gas it is simply re- 
lated to the charge q and the density n of the particles: 
A"oi = • Alternatively, for evaluating the contribu- 

tion to flux noise of pairs of particles subject to bound 
motion one may relate k to the mean distance of such 
pairs, resulting from the distribution (9): k 2 o2 — -pzy- 

This corresponds to the small-fc form of S(k) for inde- 
pendent pairs. Choosing k D \ yields a flux noise spectrum 
(figure 2) which behaves as \/oj a ( T \ where a{T) = 1 
at Tbkt and increases above. The curves for different 
temperatures cross each other (there is a hint that this 
happens in the data of reference [12]). On the other hand 
(figure 3) this temperature dependence of the exponent 
(coming from the T-dependence of the integrand in (8)) 
is compensated by using the T-dependent values of fc o2 in 
(5), see figure 3. This approach emphasizes the fact that 
large pairs, for which f(r) and l/r(r) are small, give the 
main contribution to (8) and it is the combined effect of 
the exponent a(T) and the structure factor that leads to 
a 1/oj flux noise, and all the curves for different tempera- 
tures fall on top of each other for lo > tu c , as it is seen in 
the experiments [11,12]. The white noise level showing 
up for uj < ui c is also strongly T-dependent. According 
to (5) and (7) it is given by 



S*(0) 



1 j(u> 



(2tt) 4 k B Tn v J 



0) k B T^ 2 



(13) 
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cur ves for — Iml using the combined mobility for pair 
and free motion. From more thickness to less for 
T/T B kt= 1.01,1.05,1.10. 



T he dominant contribution comes from 



' f(O) and thus 

S<f,(uj) ~ d\. However, when d c reaches the sample size L, 
7(0) saturates and the T-dependence of S<p (0) comes from 
the remaining factors in (13), in particular the density n 
of VA-excitations. This <Cmasking^> of the true critical 
slowing down has been observed in Ref. [12]. 

As a final step we combine the two contributions to 
obtain the total mobility /i(w): fi(uj) = (1 — i^)/ifc(cj) + 
VfHf(uj), with Vf = („ a ii)(d / a yi ■ Bound motion is still 



dominant up to T /Tbkt = 1.1, since its weight is larger, 
although the free mobility itself exceeds the one of bound 
pairs. Combining the two contributions (figure 4) yields a 
rather extended flat region of Im(l/e(cj)) for intermediate 
frequencies, which is another signature of <C anomalous 3> 
vortex dynamics found also in numerical simulations [24]. 

In conclusion our approximate analytical calculations, 
separating <Cfree3> and <§Cpaired;^> motion, give the fol- 
lowing insight into vortex dynamics in J J As above the 
Berezinski-Kosterlitz-Thouless temperature : 

• The anomalous Minnhagen phenomenology (MP) is 
a consquence of the motion of (unbound) vortices in a 
Coulomb potential which is screened by other <^free3> 
particles. The more this screening is reduced when ap- 
proaching the transition temperature, the more MP is 
pronounced. This regime does not lead to 1/uj flux noise. 

• V and A moving - at short enough distances and up 
to some finite life time - as a pair yield an even more 
anomalous vortex dielectric constant with temperature 
dependent frequency exponents. This effect, combined 
with a temperature dependent pair structure factor 
produces l/ui flux noise in an intermediate frequency 
range. 
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